fig_s3 = dta_overtime |> 
  mutate(accuracy_rev = 6 - accuracy) |> 
  ggplot(aes(x = date, y = accuracy_rev, color = pid, fill = pid, weight = weight)) +
  geom_smooth(show.legend = F, method = 'loess') +
  facet_wrap(~ state, ncol = 9) +
  scale_color_manual(values = c("dodgerblue3", "grey50", "firebrick3")) +
  scale_fill_manual(values = c("dodgerblue3", "grey50", "firebrick3")) +
  scale_y_continuous(breaks = 1:5, labels = c("Strongly disagree", "Disagree",
                                              "Neither agree\nnor disagree",
                                              "Agree", "Strongly agree"),
                     limits = c(1,5)) +
  geom_hline(yintercept = 3, lty = "dotted") + 
  labs(x = NULL, y = NULL,
       caption = "Missing CI segments imply CIs overlap bounds of DV scale") +
  theme_prl() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

ggsave(here("Figures", "figure_s3.pdf"), fig_s3,
       dpi = 600, units = 'in', height = 10, width = 15)
